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The spectral and thermodynamic response of segmented quantum spin chains is analyzed using a 
combination of numerical techniques and finite-size scaling arguments. Various distributions of seg- 
ment lengths are considered, including the two extreme cases of quenched and annealed averages. As 
, the impurity concentration is increased, it is found that (i) the integrated spectral weight is rapidly 

reduced, (ii) a pseudo-gap feature opens up at small frequencies, and (iii) at larger frequencies a 
discrete peak structure emerges, dominated by the contributions of the smallest cluster segments. 
The corresponding low-temperature thermodynamic response has a divergent contribution due to 
the odd-site clusters and a sub-dominant exponentially activated component due to the even-site 
segments whose finite-size gap is responsible for the spectral weight suppression at small frequen- 
cies. Based on simple scaling arguments, approximate low-temperature expressions are derived for 
the uniform susceptibility and the heat capacity. These are shown to be in good agreement with 
^ , numerical solutions of the Bethe ansatz equations for ensembles of open-end chains. 
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I. INTRODUCTION 



f| . Low-dimensional electron systems are known to be particularly sensitive to disorderEI. It is therefore difficult 
to realize pure one- or two-dimensional behavior at very low temperatures and small frequencies in nature.0 This 
is somewhat disappointing in light of recent detailed theoretical predictions for the low-enecgM scaling behavior of 
paradigmatic quantum-many-body models, such as antiferromagnetic spin chains and laddersoQ. It is the enhanced 
quantum fluctuations in these compounds, which give rise to particular low-temperature scaling regimes in very pure 
samples, and at the same instance make them highly unstable towards externally induced low-temperature transitions, 
' such as localization by impurity scattering_or three-dimensional long-range ordering due to residual small couplings 
between the lower-dimensional subsystems}!!. 

It has recently been demonstrated that in certain ID subsystems, such as random-exchange and random-spin 
chains, anomalously extended states can persist against disorderEj Eil. The physical picture is that while most spins 
are bound in randomly distributed valence bands, the unbound spins interact via virtual excitations, resulting in a 
zero-frequency band with power-law scaling.E£l Furthermore, in the case of spin ladders and spin-Peierls compounds, 
doping with randomly placed non-magnetic impurities may actually quasi-long-range ordering due to effective 

inverse-power-law interactions between the residual "pruned" spinst3't3. This replacement of an originally short- 
ranged spin liquid state by impurity-induced quasi-long-range order can be viewed as a realization of the "order by 
1 disorder" phenomenon. 

In other compounds, impurity scattering may completely destroy the connectivity within one-dimensional subsys- 
tems. If this is the case extended states cannot survive. Let us examine two specific realizations of such segmented 
, quantum spin systems: (i) CuO chains, with non-magnetic impurities, and (ii) pinned charge density waves in quasi- 
1D materials. The first situation can be realized by doping a quasi-lD compound such as SrCuC>2 with Zn. In the 
pure material, antiferromagnetic superexchange between neighboring Cu 2+ d x 2_ y 2 electrons is mediated by the filled 
2 ~ p-orbitals. By substituting Zn 2+ for Cu 2+ , static jaeancies are created, and the infinite chain is separated into 



segments of length I which follow a discrete distribution, second physical way of realizing segmented spin chains 



is the pinning of one-dimensional charge density waves. If there is competition between poorly screened long-ranged 
repulsive Coulomb forces and shortjxange attractive forces, highly inhomogeneous density wave modulations occur, 
favoring particular segment lengths .E2I Even small impurity scattering leads to a pinning of such structures. 

In both cases, there are ensembles of correlated spin segments which are most straightforwardly modelled by taking 
appropriate averages over distributions of finite clusters with open boundary conditions. The specific form of the 
distribution function strongly depends on the details of how the segments, are formed. For example, in the case of 
randomly doped CuO chains a (discrete) Poisson distribution is naturalc3u3, whereas for pinned, spatially modulated 
density waves, only two or three cluster sizes may dominate. An important factor, determining the proper distribution 
function, is the (meta-)stability of the random realizations: are they obtained from a quenched or an annealed cooling 
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procedure? While in quenched realizations all clusters that occur at high temperatures also have non-vanishing 
weight in the zero-temperature distribution function, slower "annealed" cooling processes can lead to preferred sizes 
and shapes. In particular, segments with an electronic closed shell configuration have more stable groundstates than 
others, and thus receive a higher weight in an annealed cooling process t3 

In this work, we systematically study such ensembles of antiferromagnetically correlated spin clusters, focusing 
on the evolution of the corresponding low-energy features in the dynamical spin excitation spectrum and on the 
uniform static-Sjusceptibility as a function of the hole concentration. In particular for small clusters, quantum effects 
are important^, which makes any theoretical appcpach to this problem challenging. We therefore attack this task 
numerically, using exact numerical diagonalizationE3 and scaling laws derived from conformal field theory to calculate 
the static and dynamical magnetic response for variable segment sizesE3. 

This paper is organized as follows. In the next section, a procedure to obtain excitation spectra for segmented 
Heisenberg chains is explained. The appropriate distribution functions are derived, and the evolution of the spectra 
with impurity concentration is discussed. In the subsequent section, the static magnetic response, i.e. the uniform 
susceptibility, and the heat capacity are calculated for various ensembles of finite chains. In the final section we 
conclude with a discussion of possible extensions and experimental consequences of the procedure outlined in this 
paper. 



II. EXCITATION SPECTRA AT T=0 




FIG. 1. Finite-size scaling of the staggered dynamical structure factor of the antiferromagnetic spin- s Heisenberg chain with 
open boundary conditions. The poles have been given a width of e = 0.1J. (a) chains with an even number of sites, (b) chains 
with an odd number of sites. 

To study the zero-temperature dynamical response of randomly segmented Heisenberg chains we examine the 
dynamical structure factor S(q,co,l) of finite chains with open boundary conditions: 
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5(q, w, I) = — I < n \ S «\° > ~ E n + Eo) 



(1) 



where / is the length of the finite chain, Z the partition function and n runs over all possible final states, whereas 
|0 > is the groundstate. Let us first concentrate on the staggered magnetization 



(2) 



Then S(tt,lu,1) is well defined for open boundary conditions and chains of an even or odd number of sites I. In the 
case of even I there is a unique singlet groundstate, whereas one has to take into account the doublet nature of the 
groundstate for odd I. Using numerical diagonalization techniques, we have obtained S(n, u>, I) with I = I, 20. The 
resulting spectra are shown in Fig. 1. We observe that segments of even length / have 1/2 major peaks at non-zero 
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frequencies, whereas those of odd length have a significant pole at u> = and (I — l)/2 additional peaks at higher 
frequencies. The higher energy peaks have a more complex structure, because some of the final states are singlets. 
Zero-frequency peaks occur only in the .add length segments, reflecting the fact that their groundstates are doublets. 

As in the case of closed finite chainaLij .equations derived from conformal field theory can be used to extr.act_the 
finite-size scaling behavior of S(tt, uj, 1)EB~c3 To leading order, the pole positions and amplitudes are given bytJ'B^I 

u i (l) = a i /l + (3 i /l[ln(l)}, (3) 
Ai(l) = en + bi In (1 + I) + c t In [1 + In (1 + I)], (4) 

where the coefficients aj, /3j, a*, 6j, and Cj can be treated as fit parameters. In Fig. 2 the pole positions and amplitudes 
of the lowest three poles are shown for clusters of up to 20 sites, along with the fits to the above scaling equations. 
In particular for the larger-size segments, these equations give an excellent fit to the numerical data. In the shorter 
segments higher-order logarithmic corrections for the amplitudes become increasingly relevant, and the quality of the 
fits deteriorates slightly in this regime. Note also, that the peak of the one-site cluster at ui = has an amplitude 
that does not follow the general trend. The amplitudes shown in Fig. 2 correspond to the dynamical spin response 
per site (Eq. (1)). In the following, we will consider averages over ensembles of finite chains, where the amplitudes of 
the individual segments enter as extensive quantities. In this case, the segment amplitudes per site shown in Fig. 2 
have to be multiplied by the segment length 1 (Eqs. 5 and 6). 
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FIG. 2. Finite-size scaling of the first three poles of the staggered dynamical structure factor of spin-| Heisenberg chains 
with open boundary conditions. The symbols show the values obtained by numerical diagonalizations, and the solid lines are 
the result of a fit to equations (3) and (4). Circles: first poles, squares: second poles, diamonds: third poles, (a) pole positions 
for chains of even length, (b) the corresponding amplitudes, (b) and (d) analogous results for chains of odd length. 



Now consider an infinite Heisenberg chain, doped with non-magnetic impurities, resulting in an ensemble of open-end 
segments of various lenghts 1. The average dynamical spin structure factor can then be calculated from 

S(n,u)=Y l lP(l)S(ir,u,l), (5) 
l 

where P(l) is an appropriately chosen distribution function. Using the pole structure of the response functions for 
the individual segments, we obtain 

S(tt, w ) = lP(l)Ai(l)S(cj - £*(!)). (6) 

( i 

In practice, the J-function in Eq. 6 is replaced by a Lorentzian of width e, which will be taken as e = 0.1J throughout 
the paper. P(l) determines the weight of each segment in the ensemble average, and extrinsic factors favoring certain 
cluster shapes over others enter through this function. If the chain segmentation occurs completely randomly, the 
corresponding distribution function is given by P(l) = p 2 (l — p) 1 , where p 6 [0,1] is the concentration of vacant 
sitesl!3E3. This distribution function is normalized by 

oo 

J2lP(l) = l-p, (7) 
l=i 



3 



and for p <§C 1 it can be approximated by P(l) « p 2 exp(— pi). We also note that since the total number of clusters 
per site is given by n c = J^i -f (0 = P0- ~ P)i the average length of the clusters is l av — J^i lP{l)/ n c = 1/p- The 
above distribution describes the case of quenched disorder in the infinite chain, i.e. the positions of the impurities are 
uncorrelated. In an annealed cooling process, even length segments are favoured over odd ones, because they have a 
lower groundstate energy. We describe this situation by a similar distribution function P a (l) — Cp 2 (l — p) <5o,i mod 2 
where C is determined by the normalization condition (J7|) . Only even length spin segments occur in a chain with this 
distribution function for the impurities. Note, that according to common terminology all of the ensemble averages 
which are discussed here are "quenched" because frozen disorder realizations are used. In this paper, the terms 
"annealed" and "quenched" only refer to the cooling procedures giving rise to different distribution functions. For 
true annealed disorder, however, the disorder variables would have to be treated dynamically. 

Here we wish to evaluate and analyze S(tt, u>) and S(tt) by using the finite-size scaling behavior of the lowest 
few poles in the spectrum. This few-mode approximation has been shown to be valid-Jor sufficiently large impurity 
concentrations, but it tends to underestimate the spectral response for very small pS3 Let us use the lowest three 
poles of each segment as shown in Fig. 2, i.e. the index i in Eq. (Q) runs from i = 1 to 3. Also, a cutoff length 
Imax = 10000 is used in the sum over I. The resulting frequency- integrated staggered dynamical structure factor 
S(ir) = J cLj S (iv , u) is then given by 

Zmax 

S(n) = J2lP(l)S(n,l), (8) 
z=i 

where S(n, I) is the frequency-integrated dynamical structure factor at wavevector tt, with the scaling form 

S(w,l) =a + 61n(l+0+cln[l + ln(l + /)]. (9) 

Comparing the values of S(n) obtained from exact numerical diagonalizations of ensembles of finite clusters with 
the result for S(tt) within the three-mode approximation, one finds that the three-mode approximation tends to 
underestimate S(tt) for small impurity concentrations. The difference in the integrated weight is due to neglecting the 
higher frequency poles that become more relevant for larger clusters and thus smaller impurity concentrations. This 
weight can be approximately restored by adding the properly normalized dynamical structure factor of an infinite 
Heisenberg chain, broadend by e to AS(w, u>) oc e/(ui 2 + e 2 )/7r, neglecting logarithmic corrections. It turns out that 
these corrections to S(ir) are only relevant for very low impurity concentrations and are neglegible for p > 0.2. The 
resulting spectra are shown for different p in Fig. 3 for the quenched and the annealed case. 
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FIG. 3. Evolution of S(tt,u) of the segmented antiferromagnetic spin-i Heisenberg chain as a function of the impurity 
concentration, (a)-(f) spectra for p = 0.1, 0.3, 0.5, 0.6, 0.7, 0.9. The solid lines represent the quenched case, and the dashed lines 
represent the annealed case. 

For both types of impurity distribution functions there are two main features that occur in S(n, u) upon increasing 
the number of vacancies in the infinite chain. First, the integrated spectral weight decreases rapidly upon increasing 
the impurity concentration. This is shown explicitly in Fig. 4, where the solid lines in part (a) and (b) represent the 
dependence of S(tt) on p. Note that in the annealed case the integrated spectral weight is decreasing at a slower rate 
than in the quenched case, especially for large impurity concentrations. This is due to the fact, that the integrated 
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spectral weight of finite chains is a monotonous function of I. Thus the total response increases if at constant impurity 
concentration the odd- length segments (starting with 1 = 1) are substituted by even-length segments (starting with 
1 = 2). Most of the suppressed spectral weight comes from the low-energy continuum of the response function due to 
the large segments. It can be shown by considering a lowest-order single-mode approximation that S(tt,oj) is in fact 
exponentially suppressed at low frequencies^. 

The second common feature is the emergence of a discrete peak structure at larger frequencies (of order J), dominated 
by the contributions of the smallest cluster segments. The dominant segments, occuring according to the distribution 
function P(l), can be identified from these higher energy peaks. In the annealed case, the smallest segments are the 
two-site clusters with a pole in S(ir,U),l) at oj = 1J, and the four-site chain with a major pole at u) = 0.66 J. These 
poles are well separated from the low-energy continuum and carry most of the spectral weight of S(tt, oj) at impurity 
concentrations p > 0.5 (Fig. 3 (c) - (e)). The pole of the smallest segment at uj = 1J dominates the dynamical 
response function in the annealed case at high concentrations, p > 0.7. One thus expects that most of the spectral 
weight can be obtained from the two-site cluster for p close to one. In fact, from Eq.(8) one finds 

5(7r)«2P(2)5(7r,2)«i^, P ->1, (10) 

explaining the linear behavior of S(tt) for large impurity concentrations. Since in the quenched case the single spin 
sites and the three-site chains are also present, their low frequency poles contribute strongly to the spectral weight. 
The major pole of the three site chain at uj = 1.5J is the first well-defined pole to separate from the low-energy 
continuum upon increasing the impurity concentration, as can be seen in Fig. 3 (a) already at p = 0.1. It dominates 
the spectrum away from the low-energy pseudo-branch until, at intermediate concentrations, the lowest-energy poles 
of the two- and four-site clusters are also separated from the continuum. Characteristic features of the underlying 
impurity distribution can thus be read off from the higher frequency part of the spectrum. 

The major difference between the two impurity distributions we have studied lies in the rate at which S(it,uj) is 
suppressed for small frequencies. Consider the annealed case first. The even-site segments do not have a pole at 
uj = due to their inherent finite spin-gap. Thus the exponential suppression of the contributions from the large 
segments leads to the development of a pseudogap at small frequencies with increasing impurity concentration. The 
finite values of S(tt, 0) are due to our replacement of the 5-peak in Eq. (6) by Lorentzians. This mimics the various 
broadening mechanisms in real materials, such as thermal broadening, scattering from phonons, and interactions with 
out-of-chain impurities. As shown in Fig. 4 (c), the residual spectral weight S(ir,0) is exponentially suppressed with 
increasing p. In Fig. 4 (d) the important higher-energy peaks are compared to S(tt, 0), clearly indicating the reduction 
of the zero-frequency weight in the annealed case. 
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FIG. 4. (a) and (c) The frequency-integrated dynamical structure factor at wavevector n, S(n) (solid line) and the residual 
spinon density of states S(n, 0) (dashed line) as a function of the impurity concentration p. (a) quenched disorder, (c) annealed 
disorder, (b), (d) Dependence of the amplitudes of the peaks at uj = (full circles) uj\ = 0.66J (solid line) , UJ2 = 1.0J (dashed 
line) and LJ3 = 1.5J (dotted line) on the impurity concentration p. (c) quenched, (d) annealed. 

Now consider the quenched case. All odd-length segments and especially the dominating one-spin segment show a 
pole in their dynamical response function at zero frequency. Therefore one expects the rate of suppression of S(tt, 0) 
to be significantly reduced in the quenched case with respect to the annealed case. This is clearly observed in Fig. 
4, comparing the value of S(ir, 0) for both distributions. In the quenched case the weight at uj = remains a major 
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contribution to S(n), even at large impurity concentrations. Fig. 4 (a) shows that the reduction of S(n,0) with 
increasing p is less pronounced in the quenched case. In fact, the amplitude of the peak at u> — dominates the 
spectrum for all p, as shown in Fig 4. (b). 



III. THERMODYNAMIC RESPONSE: UNIFORM SUSCEPTIBILITY AND HEAT CAPACITY 



In this section we examine the thermodynamic response of segmented spin-1/2 Heisenberg chains. Based on simple 
scaling arguments, low-temperature approximations for the uniform susceptibility and the specific heat are derived 
and shown to be in good agreement with a numerical solution of the problem, using the Bethe ansatz equations for 
open-end chains. The impurity tail of the low-temperature susceptibility (i.e. the divergence of x(T) at T = 0) in an 
ensemble of segmented chains is caused by the odd-length segments. Before examining the sub-dominant contributions 
of the even-length segments it is therefore necessary to first identify and discuss the dominant divergent contributions 
of the odd-length clusters. This can be achieved by analyzing [x(T)T]. Its value at zero temperature, \imT^o[x(T)T], 
gives the prefactor of the low-temperature impurity tail in x(T)- For chains of even length, which always have a finite 
spin gap, the value of [x{T)T] approaches as T goes to zero, whereas for odd lengths i, \m\T->o[x{T)T] = 1/(4/). It 
follows that in the quenched case 

[x(T)THp) = lim[ X (T)T](p) = ldmf>(0[x(T)T](i) = (H) 

with a maximum at p = 2 — \/2 w 0.586. In the ideal annealed case only even-length segments contribute to the 
average, and therefore [x{T)T] (p) vanishes. The Curie constant is C(p) = lim' r _ 00 [x(r)T] = (1 — p)/4, independent 
of the length of the contributing segments. 
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FIG. 5. Linearized form of the scaling behaviour for the low temperature susceptibility of segmented Heisenberg chains for 
impurity concentrations p = 0.05 (top) ,0.1,0.2,...,0.9 (bottom). 



To study the temperature dependence of x(T) we can calculate the finite-size susceptibilities from the Bethe ansatz 
equations which determine the energy levels of open Heisenberg chains: 



Xk + 1 
Xk - i 



M 

n 

3+k 



Xk — Xj + 2i Xk + xj + 2i 



(12) 



Here the number of roots M determines the total S z component of the state through the relation S z = 1/2 — M. 
In logarithmic form this equation becomes 



M 



21 ctan x (xk) = + /J 



ctan 



— l / *^ 



ctan 



1 / Xk T" Xj 



(13) 



G 



where all the Ik are integers with k = 1, M. Given a solution of the above equations for a set of I^s, the energy 
of the corresponding eigenstate is 



E 
1 



1 1 M 
— -2V 
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fc=l 



1 



(14) 



The groundstate in a given S 12 sector, Eq(S z , I), is obtained from the set {Ik} = {l+l, 1+3, 1+5, 1+2AI — 1}. And the 
first exited state in the S z > sectors can be determined from the set {I k } = {2S Z , 2S Z + 3, 2S Z + 5, 2S Z + 7, ...,^-1}. 
Using an iterative method to determine the energy gaps for the even length segments, AEi — Eq(S z = 1,1) — Eq(S z = 
0, 1), one obtains for the low temperature susceptibility of the segmented chain: 



[x(T)T](p) = [x(T)TUp) + 2^2' T 



P(l) 



1=1 



exp(A£7j/ r )' 



(15) 



where the sum is restricted to even I. This equation can be evaluated numerically, using the energies from the finite- 
size Bethe ansatz equations (Eq. 14) and an appropriate distribution function P(l). It is found that the results for 
[x(T)T] typically converge below a length cut-off around l m ax — 1000. Results for the quenched case are shown in 
Fig. 5. For low temperatures there appears to be a universal scaling behavior which can be elucidated by expanding 
Eq. 15. 

To determine the low-temperature scaling behavior of the susceptibility we consider the lowest-order finite-size 
scaling behavior of the energy gap, AEi w a±/l, neglecting higher-order logarithmic corrections (Eq. ||) which 
become increasingly important when larger segments occur (i.e. at very small impurity concentrations). For low 
temperatures one obtains 



[x(T)T](p) - [x(T)T] (p) a / exp 

lo 



dl 



4ai 



Ki 



I 4ai . / 1 
■ In 




In- 



exp 



The expected low-temperature behavior is thus of the form 
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(16) 



where 7 turns out to be 



[\a)T][^-[ X iTa- A P) x (j\ > 




(17) 
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FIG. 6. Fit of the effective low-temperature gap 7. The circles show the numerical solution of the Bethe ansatz equations, 
and the solid line is the best fit with on = 3.88J. 
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The observed linearity in Fig. 5 where we have plotted \n[( X (T)T(p) - [ X (T)T} (p))(T/ J) 1 / 4 ] vs. (T/J)- 1 ' 2 for 
various p confirms the validity of this low-temperature expansion. The free parameter a\ can be determined by fitting 
j(p) of Eq. |l8|to the exact numerical Bethe ansatz solution of [x(T)T](p] (Eq. 15). As shown in Fig. 6, for sufficiently 
small concentrations the fit turns out to be very good, giving ot\ — 3.88 J. 

A similar procedure can be used to determine the scaling behavior of the specific heat of segmented Heisenberg 
spin-1/2 chains. For the quenched case, one finds 

Coc (j\ VV/V5VJ (19) 
at sufficiently low temperatures. The corresponding effective gap 7' turns out to be 



7' = 70 + ^/4 ^(ji-). (20) 
where a fit of the numerical data yields 70 = 0.182 and a[ = 4. 106 J. 



IV. CONCLUSIONS 



In summary, we have examined the spectral and thermodynamic response of various ensembles of segmented an- 
tiferromagnetic Heisenberg spin-1/2 chains. We have calculated the dynamical spin structure factor which can be 
probed by Neutron scattering experiments. The particular dependence of this quantity on the impurity concentration 
is determined by the distribution of segment lengths. However, there are several generic features which are observed 
for the most common distribution functions. These are (i) a rapid decay of the integrated spectral weight with in- 
creasing impurity concentration, (ii) a suppression of low- frequency poles (pseudo-gap), and (iii) the emergence of a 
discrete pole structure at higher energies, dominated by the smallest contributing clusters in the average. While we 
find that the first two predictions are consistent with the presently available experimental data, point (iii) might be 
the hardest to verify because the corresponding signals in a highly disordered sample are typically rather small and 
broad. 

Two main contributions to the low-temperature thermodynamic response of segmented chains can be identified: a 
dominant divergent component from the odd-side clusters and a sub-dominant exponentially activated component due 
to the even-site clusters. The sub-dominant contribution is analyzed by-sutttracting the impurity tail from the total 
response function, as it is commonly done in the analysis of expcriments.oEEl Assuming a certain type of distribution 
function, the corresponding effective gaps, here 7 for the uniform susceptibility and 7' for the specific heat, can be 
calculated. An analysis of these gaps can thus be used as an indicator of the underlying distribution function of 
segment lengths. 

Complete segmentation of quantum spin chains, as it has been treated in this work, can be viewed as an extreme 
case of impurity scattering with an infinitely large on-site repulsive potential. While this mechanism may indeed 
lead to segmentation, in many physical realizations longer-range exchange paths exist which can partially restore 
extended states of the undoped parent systems. Consider for example Zn-doped CuGe03. This compound is known 
to have sizeable next-nearest neighbor exchange interactions along the Cu02-chain direction, giving rise to an effective 
Ji — Ji model. Below the transition temperature T$p the compound goes into a dimerized spin-Peierls phase, whereas 
above T$p it is in a critical quasi-one-dimensional state. This spin-Peierls transition is suppressed upon replacing 
the Cu-atoms randomly with ZrO, most likely because the partial segmentation due to the non-magnetic impurities 
impedes the quantum-critical extended states within the chains which in turn are a prerequisite for Peierls transitions. 
Because of the longer-range exchange paths J2, these extended states are not completely destroyed, and a remnant 
spin-Peierls phase is observed in Cui_ x Zn x Ge03 at sufficiently small impurity concentrations x.Ej 

The segmentation of critical one-dimensional system competes with transitions such as three-dimensional order- 
ing due to small inter-chain interactions or Peierls- type transitions. Segmented one-dimensional phases are more 
stable against low-temperature ordering transitions, and critical states can in turn be created in a controlled man- 
ner by impurity-doping short-range ordered systems, for example by. intrpd ucing non-magnetic sites into quasi-one- 
dimensional spin liquids such as two-leg spin-1/2 Heisenberg ladders .llSEJcJ 

We thank A. Honecker and B. Normand for useful discussions, and acknowledge the Zumberge Foundation for 
financial support. 
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